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Abstract 

It has been known through some examples that parameters of an electromagnetic 
medium can be so defined that there is no dispersion equation (Presnel equation) to 
(y-j : restrict the choice of the wave vector of a plane wave in such a medium, i.e., that the 

dispersion equation is satisfied identically for any wave vector. In the present paper, 
a more systematic study to define classes of media with no dispersion equation is 
attempted. The analysis makes use of coordinate-free four-dimensional formalism 
in terms of multivectors, multiforms and dyadics. 
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1 Introduction 

Time-harmonic plane waves in linear (bi-anisotropic) electromagnetic media are charac- 
terized by dispersion equations [1] (or Fresnel equations [2]) restricting the choice of the 
k vector of the plane wave for any given frequency u. In the general case the dispersion 
equation is of the fourth order. For a restricted class of media the fourth-order equation 
can be reduced to two second-order equations, in which case the field solutions can be 
decomposed in two independent sets obeying certain polarization properties [3, 4, 5, 6]. 
A simple example of such a medium is the one defined by uniaxial permittivity and per- 
meability dyadics. In a still more restricted case, e.g., that of isotropic media, the two 
equations coincide to a single second-order equation, whence the medium appears free of 
birefringence [7, 8, 9, 10]. Finally, one can define media for which the dispersion equation 
is satisfied identically. In such a case, the medium imposes no restriction on the choice of 
the k vector of the plane wave. It appears that this property has been widely neglected in 
the literature. However, media belonging to this cathegory have recently emerged in vari- 
ous studies [11, 12, 13, 14]. Interestingly, it has been shown that novel types of boundary 
conditions may arise at the interface of some media of this kind [12, 15, 16]. It is the 
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purpose of the present paper to define, in a more systematic manner, classes of media in 
which fields are not restricted by a dispersion equation. 

Because of its compactness, the analysis will apply the four-dimensional differential- form 
formalism of reference [17] empowered by coordinate-free dyadic representations for linear 
mappings [18]. The source- free Maxwell equations can be expressed in simple form as 

dAf = 0, (1) 
dA* = 0, (2) 

where the electromagnetic two-forms $, f G F 2 are defined by 

$ = B + EAe 4 , (3) 
* = D-HAe 4 , (4) 

in terms of the spatial 3D one-forms E, H e F x and two- forms B, D. The temporal 
one-form 

£4 = dr, r = ct, (5) 

can be completed by spatial one-forms £±,£2, £3 to a basis of one-forms. The comple- 
mentary basis vectors e 1 ,e 2 ,e 3 ,e 4 e Ei satisfy e^Sj = 5ij. Details of the notation and 
operational rules applied in this paper can be found in [18]. For similar equations in 
tensorial notation one may visit the references [2, 14]. 



2 Dispersion dyadic 

A plane wave can be defined by fields depending on the space-time vector x = r + e 4 r as 

*(x) = */Mx), (6) 
*(x) = tf/Mx). (7) 

Here, /(.) is any differentiable scalar function and v G Fi is the wave one-form which 
contains the wave vector k and the frequency uj of the Gibbsian 3D vector notation. For 
a plane wave the Maxwell equations (1) and (2) imply 

= 0, (8) 

j/A* = 0. (9) 

Applying the identity [18] 

a J(i/ A*) = v A (aj*) + (a|i/)*, (10) 

and similarly for ^, valid for any vector a, assuming a\v 7^ 0, we can express the field 
two-forms in terms of potential one-forms <f>, ip as 

$ = uA <f>, 4>=-A-, (11) 

a\v 

a|* 

* = i/A-0, ip = j —. (12) 

a \u 
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A linear medium can be represented in terms of a linear mapping between the electro- 
magnetic two- forms. Let us assume that the medium is homogeneous and time- invariant, 
i.e., independent of x. In terms of a bidyadic M e F 2 E 2 , mapping two-forms to two- forms, 
one writes _ 

* = M|$. (13) 

Another, equally valid representation, makes use of another bidyadic N e F 2 E 2 as 

$ = M|*. (14) 

M and N may be expanded in any given bivector basis and two-form basis as 6 x 6 matrices, 
each involving 36 scalar parameters in the general case [19]. When the two bidyadics are 
of full rank, they are inverses of each other, while in the converse case only one of the 
two representations (13) and (14) may exist. The medium equations can be equally well 
expressed in terms of modified medium bidyadics M m , N m e E 2 E 2 as bivector equations 

e^L* = M m |$, (15) 
e N [Q = N m |tf, (16) 

in terms of the quadrivector = ei 234 and 

M m = e N [M, I m = e^LN (17) 

mapping two-forms to bivectors. 

According to Hehl and Obukhov [2], the general medium bidyadic can be decomposed in 
three components as 

M = Mi + M 2 + M 3 , (18) 

respectively called the principal, skewon and axion parts of M. The axion part is a multiple 
of the unit bidyadic, M 3 = M 3 T( 2 ) T , while Mi and M 2 are trace free. The skewon part is 
characterized by the modified bidyadic M m2 = ejv[M 2 which is the antisymmetric part of 
M m , while the principal modified bidyadic M ml = e^Mi is symmetric. The other medium 
bidyadic N can be given a similar decomposition. 

Condition (9) defines an equation for the potential one-form <fi as 

i/A* = (i/AM|_i/)|0 = O. (19) 

The dyadic in brackets belongs to the space F 3 E X mapping one-forms to three-forms. An 
equivalent equation can be formed in terms of the dyadic D(u) G EjEi mapping one-forms 
to vectors _ _ _ _ 

D(i/) = e N [(iy AM[u) = -u\ (e N [M) \y = M m [\uu, (20) 

and called the dispersion dyadic, as 

B(v)\<f> = 0. (21) 
From the form of (20) it follows that the dispersion dyadic also satisfies 

T5(v)\v = (22) 
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for any v. Assuming non- vanishing field two-forms <fr = v A <fi 7^ 0, the one- forms <fi and 
v are linearly independent, whence the rank of the dispersion dyadic is demanded to be 
less than 3. Thus, D(u) must satisfy the equation [18] 

PV) = ^{u)fi5(u)fi5(v) = 0, (23) 
6 

for any v characterizing a plane wave in such a medium. One can show that we can 
expand 

W\v) = (e N [u)(e N [u)D(u) = 0, (24) 

so that (23) is actually equivalent to a scalar equation which can be given the following 
explicit form [20], 

D{u) = y N s N \\(M m 2 \v)) 

= ^e N s N \\(M m ^iy\\(%^(iyiy\\U m )))) = 0. (25) 

(25) (or its equivalent dyadic form (23)) is called the dispersion equation, because for 
time-harmonic fields it defines the relation between the wave number and the frequency 
of the wave. Alternatively, it is called the Fresnel equation. Obviously, (25) is of the 
fourth order in v. Other forms for the dispersion equation are found in [21, 22, 2, 23, 14] 

Starting from the medium condition (14), a similar equation in terms of the bidyadic N m 
is obtained, 

D\v) = ^^||(M^(^JJ(M m ^(^JJ>1 m )))) = 0. (26) 

When the medium bidyadics M, N are of full rank, the ensuing two dispersion equations 
must have the same solutions since they correspond to the same plane wave [10]. In 
particular, if one of them is satisfied identically for all u, so must be the other one. When 
the rank of M is less than six, the bidyadic N and the equation (26) do not exist, and 
conversely. 



3 Media with no dispersion equation 

To define media with no dispersion equation, either (23) or (25) must be identically valid 
for any one-form v. Two examples have emerged in the past, the skewon-axion media 
(IB-media) [2, 11] and the P-media [13, 14]. Actually these two media will also come out 
through the present analysis, added by a third one denoted by case 1 below. 

Starting from (23), requiring that the dispersion dyadic be of rank less than three, it 
follows that there must exist four vectors in terms of which we can write 

D(i/) = ac + bd. (27) 

This corresponds to the case of planar dyadics in the 3D Gibbsian formalism [24, 25]. 

Since the left side of (27) is a quadratic function of u, so must be the right side. This 
gives us a few possibilities. 
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1. Each of the four vectors is a linear function of v 

2. a and b are quadratic functions while c and d do not depend on v 

3. a and d are quadratic functions while b and c do not depend on v 

4. a is a quadratic function while b and d are linear functions and c does not depend 
on v 

Other possibilities do not seem to bring any new solutions because the dispersion equation 
(25) is invariant to replacing the modified medium bidyadic M m by its transpose M^. In 
fact, this implies replacing the dispersion dyadic D(i/) by its transpose D T (u) and D^(u) 
by D^ T (u). Let us consider these four cases separately. 

3.1 Casel: a, b, c and d are linear functions of v 

Because the dispersion dyadic satisfies D(i/)|i> = u\D(u) = for all u, the four vectors 
must satisfy 

a(c\u) + b(d\u) = 0, (i/|a)c + (i/|b)d = 0. (28) 

Assuming 

U (2) = (aAb)(c Ad) ^ 0, (29) 

which requires that the dispersion dyadic be of rank 2, the vectors a, b must be linearly 
independent and so must c, d, whence we obtain 

&\v = h\u = c\u = d\v = 0, (30) 

for all v. This requires that there must exist bivectors A, B, C, D such that we can 
express 

a = A[^, b = B[^, c=C|_i/, d = D[^, (31) 

and _ _ 

D(i/) = M m [ \yv = (AC + BD) |_ \yv. (32) 

Because 

(M m -(AC + BD))LL^ = (33) 

must be valid for all u, applying the result given in Appendix 1, the modified medium 
bidyadic M m must be of the form 

M m = AC + BD + ae N [T (2)T , (34) 

for some scalar a. Equivalently, the medium bidyadic M must be of the form 

M = IIC + AD + aP )T , (35) 

with n = e N [A and A = e N [B. The expression (35) defines one possible class of media 
in terms of two bivectors, two two-forms and a scalar, for which the dispersion equation 
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is satisfied identically. The medium bidyadic N = M 1 has a similar form, which can be 
shown by forming the inverse bidyadic of (35), details of which are given in Appendix 2. 

In the special case D^(i/) = for all u, (34) is reduced to the form 

M m = AC + ae N [\^ T , (36) 

while requiring D(u) = leads to M = M~\^ T , i.e., a medium with just an axion com- 
ponent. One may note that the medium defined by (34) includes any Q-medium [18], 
defined by 

M m = MW\ (37) 

when the dyadic Q e EiEi is antisymmetric. In fact, because such a dyadic can be 
expressed as Q = A[T T , where A is some bivector, expanding 

M m = M(ALF) (2) 

= M(AA-^v|(AAA)T( 2 ) T ), (38) 

yields a special case of (36). 

One may note that the condition (29) requires that none of the four bivectors in (34) is 
simple. In fact if, for example, we assume that A simple, whence we can express it in 
terms of two vectors as A = ai A a 2 , we have a = A\y = (a 2 |f)ai — (ai|i/)a 2 = for any 
v satisfying &y\v = and a 2 |f = 0. Thus, (29) is not valid for one-forms v taken from 
the corresponding 2D subspace. In such a case the dispersion dyadic D(u) is of rank 1 for 
a 2D subspace of v and of rank 2 for other one-forms v. Since media defined by (34) have 
no dispersion equation regardless of one or more of the bivectors being simple, dispersion 
dyadics D(v) of variable rank are included in addition to those of rank 2 or rank 1 for all 
v. One may conjecture that (34) represents all solutions for case 1 media. 

3.2 Case 2: a and b are quadratic, while c and d are independent 
of v 

Let us now assume that in (27) the vectors c and d are independent of v, whence a and 
b must be quadratic functions of v. In this case it appears preferable to start from the 
representation 

v A = -e N [D(is) = Tc + Ad, (39) 

where the two three-forms 

r = -e N [sL, A = -e N [b, GF 3 (40) 

are now quadratic functions of v. Since any two linearly independent three- forms can be 
interpreted as belonging to a basis of three-forms, which can be constructed in terms of 
a basis of one-forms, say a, (3, 7, 6, we can express 



r = aA/3A7, A = a A/3 AS. 



(41) 
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Thus, (39) can be written as 

i/AM|_i/ = aA/3A(7C + £d), (42) 

or, more generally, as 

v A M[u = a A /3 A P T . (43) 

Since terms of the form aa' + /3b' can be added to the dyadic in brackets in (42), we 

may well assume that the dyadic P G EjFi in (43) is of full rank. The special case 

corresponding to D^ 2 \u) = 0, i.e., when the bracketed dyadic is reduced by assuming 
dd = 0, will be considered separately. 

Before proceeding, let us briefly check that when M satisfies a relation of the form (43), 
the dispersion equation is satisfied identically. The dispersion dyadic can be expressed in 

the form 

D( u ) = e N [(y A M [u) = D [P T = (D J T ) | W, (44) 

where the bivector D G E 2 is defined by 

B = e N [(aAf3), (45) 

and D [l T G EjEi is an antisymmetric dyadic. The bivector D is simple because it satisfies 

£ W |(DAD) = e N \((e N [(a A P)) A (e N [(a A (3))) = e N \(a A (3 A a A (3) = 0. (46) 

Since a simple bivector can be expressed in terms of two vectors as D = d x A d 2 , the 
antisymmetric dyadic satisfies 

(DLPY 3 ) = ((d 1 Ad 2 )LF)( 3 ) = (d 2 d 1 -d 1 d 2 )( 3 ) 

= i(d 2 di - d 1 d 2 )^(d 1 d 2 ^d 2 d 1 ) = 0. (47) 

From this and the rule (A|B)( p ) = A( p )|B( p ) it follows that the dispersion equation 

B( 3 )(i/) = (DLF) (3) |P (3)T = (48) 

is, indeed, satisfied for any v and any medium bidyadic M satisfying (43). 

Returning to (43), two conditions for the quantities on its right-hand side are obtained 
by noting that the left-hand side vanishes when multiplying by \v from the right or by 
vA from the left. These lead to the respective conditions 

ol A (3 A (F^i/) = 0, (49) 

and 

i/AaA/3AP T = 0. (50) 

Condition (49) requires that the three one-forms ct, f3 and P T \v be linearly dependent, 
whence they must satisfy a relation of the form 

Aa + Bp + C(F T \i/) = 0. (51) 
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For C = we have a A (3 = and D(u) = corresponding to the axion medium. For 
C 7^ we can choose C = — 1, whence the relation can be rewritten as 

f T \v = Aa + B(3. (52) 

Omitting vanishing of both A and B which would again yield the axion medium, without 
sacrifying generality (a and (3 appear interchangeable at this stage), we may set B — 1, 
whence 

/3=P T \is- Act. (53) 
Inserting this in (43), the condition becomes 

v aM|_i/ = a A (F» AW. (54) 

Since the right-hand side of (54) must be a quadratic function of v, it appears that P 
must be independent of v. Thus, ex must be a linear function of v and it can be expressed 
as 

a=W\v, (55) 
in terms of some dyadic B e E]Fi. Inserting this in (54) yields the representation 

uaW[u= (W\v) A(W\v) AF T , (56) 

Applying now the second condition (50), i.e., multiplying (56) by vA from the left, leads 
to the condition 

v A (B T |i/) A (P T |i/) A P T = 0, (57) 

valid for any one-form v. Recalling that P was assumed to be a dyadic of full rank, we 
can multiply (57) by |P _1T from the right. The condition is then reduced to 

v A (Tf T |^) A (F T |i/) A F = u A (W~\v) A (F T |i/) = 0. (58) 

Since the three one-forms are linearly dependent for all v, the dyadics must be related as 

A] + B^ + CP = 0, (59) 

for some scalars A, B, C. This leaves us the following three possibilities: 

1. A = which implies BB T \v + CP T \v = 0. Since v A M \y = 0, we must have 

M = a\^ T . (60) 

2. A and B = implies F = -(A/C) = \ and v aM\v = (A/C)u A (Efy T \y. In this 
case we must have _ _ 

M = (B£T) T + aP^ T , (61) 

3. A ^ and B ^ implies B = -(A/B) = \ - (C/B)P. In this case we must have 

uaW[v= -{A/B)u A F t |i/ AF = -{A/B)u A W )T [u, (62) 

and 

M = P^ T + aP^ T . (63) 
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In (61) and (63) the coefficients A,B,C have been suppressed. Thus, in case 2 the media 
with no dispersion equation fall into two classes, those of skewon-axion media [11] and 
P-axion media while axion media are special cases of both of the two. The class of P-axion 
media is a generalization of that of P-media [13], with an added axion term. 

Actually, all of the above medium bidyadics can be expressed under the single, more 
general, form as 

M = + B^\) T + Cl^ T (64) 

where A, B, C are any scalar coefficients. In fact, (64) yields the skewon-axion medium 
for A = and the P-axion medium for A ^ 0. For A ^ this can be seen by defining 
new symbols as 

P = B + -!, M — A, a = C -—. (65) 

Let us finally test whether (64) really satisfies the dispersion equation (23) for any one- 
form v by expanding the dispersion dyadic as 

D(i/) = -e N \{y AW\y) 

= -e N [{Av aW^[v + Bv A$£\) T \y) 

= -e N [{Au A (W\u) aW + Bu A A F) 

= (-e N [(uA(W\iy)))[(A§ + B]) T 

= (A[F)|(Al + BT) T . (66) 

Here we have denoted A = — e N [(u A (B T |z/)), which is a simple bivector satisfying 
A A A = 0. Now the dispersion equation (23) becomes 

W 3 \u) = (ALH (3) |(Al+ B = \)^ t = 0, (67) 

which is valid for the same reason as (47). It is remarkable that this is valid for any dyadic 
B which means that the assumption of full-rank dyadic P can be relaxed. 

The previous analysis assumed D^(u) ^ for all v. The converse case can be handled 
by replacing (27) by 

= ac, (68) 
where a is a quadratic function of v and c is independent of v. From 

U(i/)|i/ = a(c|i/) =0, (69) 

satisfied for all v, we must have either a = or c = 0, both of which lead to D(u) = 0. 
This, again, corresponds to the axion medium case which is included the solution (64). 



3.3 Case 3: a and d are quadratic in v while b and c are inde- 
pendent of v 

In this case we assume that in the dispersion dyadic representation (27) the vectors c and 
b are independent of v. From 

b A 15{y)\v = b A (ac + bd)|i/ = (bAa)(c|i/) = 0, (70) 
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valid for all v, we conclude that either c = or b A a = 0, whence a must be of the 
form a = (A ||ff)b for some dyadic A e EiEi. Similarly, we must have either b = or 
d = (D ||f^)c. Thus, the dispersion dyadic of case 3 coincides with that of case 2, or its 
transpose, with D^ 2 \u) = for all v. 

3.4 Case 4: a is quadratic, b and d are linear and c is indepen- 
dent of v 

In this final case from 

U(i/)|i/ = a(c|i/) + b(d|i/) = 0, (71) 

assuming aAb ^ we obtain c\v = and d\v = for all v. This requires c = and 
d = D \y, whence the dispersion dyadic must be of the form 

U(j/) = BDLL^ (72) 

for some bivectors B, D This leads to case 1 solution (34). Alternatively, from a A b = 
we obtain the relation a = Ab = (a'|i/)b, since A must be a linear function of v. In this 
case the dispersion dyadic becomes 

=b((a»c + d), (73) 

which reduces to the form (72). In conclusion, case 4 coincides with the special case 1 
satisfying D^fa) = for all v. Thus, the corresponding modified medium bidyadic has 
the form of (36). 



4 Discussion 

The previous analysis shows us that the possible medium bidyadic M, in which a plane 
wave is not rectricted by a dispersion equation, is either of the general form (34) or (64). 
Starting from the representation (14), similar expressions can be obtained for the bidyadic 

N, which equals M^ 1 when the two bidyadics are of full rank. It is shown in Appendix 
2 that both of these representations lead to the same class of media for case 1. Let us 
now study whether the same property is valid for case 2. It is known from previous work 
that for pure axion, skewon and P-medium bidyadics the inverse bidyadics, when they 
exist, have the same axion, skewon or P-medium character [11, 13, 14]. So, one may 
wonder whether this property will be valid for skewon-axion and P-axion media as well. 
Since there are no analytic expressions known to these authors for the inverse of general 
skewon-axion or P-axion bidyadics, the problem must be considered through some known 
properties. 

It is known that the modified medium bidyadic of a P-medium 

M m = M eiV p 2)T , (74) 
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for some dyadic P G E1F1 satisfies the quadratic equation [26, 5, 27] 

M^-M m = p ejv LT (2)T , P = M 2 tiW\ (75) 

The natural dot product between two bidyadics A, B or two two-forms <J?, \I/ is defined 
by 

A B = A|(£jvJ (2) |B = B- A, (76) 
= $|( eiV J (2)T |* = * • *• (77) 

The converse is not necessarily true, since, as was shown in [5], a solution of (75) can 
alternatively be a bidyadic of a Q-medium, which has the form [18] 

M m = MW 2 \ (78) 

for some dyadic Q G EiEi. However, as is shown in Appendix 3, if a solution of (75) is of 
full rank (P ^ 0) and belongs to case 2, it must be a P-medium bidyadic, since Q-medium 
bidyadics possessing no dispersion equation fall necessarily into case 1. 

From (75) it follows that the bidyadic of a P-axion medium 

M m = Me^^ (2)T + «e w J (2)T , (79) 
satisfies the bidyadic second-order equation 

• M m - a(W m + M m ) = (P- a 2 )e N \} {2)T '. (80) 

Let us multiply (80) by M m 1T | from the left and by (M" 1 from the right, whence we arrive 

at 

s N C\ {2) -«(M m 1T + M" 1 ) = (P-« 2 )M- 1T |(e JV Ll (2)T )|M- 1 . (81) 

For P 7^ a 2 we can find the equation corresponding (80) for the bidyadic N m = ejvejvLLMm 1 ' 

K-%n+ p^rc + ^n) = J (2)T (82) 

Since (82) is of the same form as (80), we may conclude that, in the general case, the 
inverse of a P-axion medium is another P-axion medium. 

In the special case when P-axion medium bidyadic (75) is restricted by P = a 2 , i.e., for 

M m = Me N [P (2)T + ae N J (2)T , a = ±VP, (83) 
with either sign, (82) reduces to 

a(C + M m ) = ejv LT (2)T , (84) 
which is satisfied by bidyadics of the form 

M m = l+^-e N \l {2)T , (85) 
la 



12 



M 


N 


axion 


axion 


skewon 


skewon 


P-medium 


P-medium 


skewon- axion 


special P-axion 


special P-axion 


skewon-axion 


general P-axion 


general P-axion 



Table 1: Relations between various case 2 medium bidyadics. 



where A is any antisymmetric bidyadic. From this we conclude that the inverse of a 
special P-axion medium (83) is a skewon-axion medium. This works also backwards so 
that the inverse of a skewon-axion medium is of the special P-axion medium form (83). 
This is, however, not valid for a — 0, which corresponds to the pure skewon medium 
whose inverse is another pure skewon medium [28], while for P = the bidyadic M m in 
(83) does not have any inverse. The result is shown in Table 1. 

It may be of interest to note that both case 1 and case 2 classes of media with no dispersion 
equation are closed in affine transformations. Let us consider a full-rank dyadic A e E 1 F 1 
mapping vectors as x a = A|x, bivectors as A a = A^|A and medium bidyadics as [18] 

M a = A(- 2 ) T |M|I( 2 ) T (86) 

Applying this transformation, case 1 medium bidyadics are mapped as 

M a = A^Knc + AD W (2)T )|I (2)T 

= II a C a + A a D a + aT {2) , (87) 

with 

n a = A(- 2 ) r |II, A a = A(- 2 ) T |A, (88) 

C a =A( 2 )|C, D a = A( 2 )|D, (89) 
and case 2 medium bidyadics are mapped as 

M a = I^ t \(aW^ t + b(M) t + cP^ t \1^ t 

= AW^ T + B{^J) T + CP^ T , (90) 

with _ 

B a = A|B|A x . (91) 

Thus, in both cases, the form of the medium bidyadic remains the same after the trans- 
formation. 



5 Conclusion 



This paper concerns the problem of defining classes of electromagnetic media in which 
plane waves are not restricted by a dispersion equation (Fresnel equation). Some media 
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with such a property have recently emerged in various studies. In the analysis four- 
dimensional formalism in terms of multivectors, multiforms and dyadics was applied as 
based on reference [18]. 

Requiring that the dispersion dyadic D(f) satisfy the condition D^ 3 \u) = for all possible 
wave one-forms v the problem was split in four cases which were analyzed separately. Case 
1 and case 2 were shown to yield medium bidyadics M, respectively defined by (35) and 
(64), with the desired property, while cases 3 and 4 could be reduced to cases 1 and 2 with 
no new solutions. Previously known media showing no dispersion equation were seen to be 
special cases of the results obtained here. In particular, the Q-medium with antisymmetric 
Q dyadic belongs to the class of case 1 media (35) while the skewon-axion medium and 
P-medium are both special cases of the class of case 2 media. When considering the 
problem in terms of the alternative medium bidyadic N, no new solutions could be found 
in any of the cases considered. 

Appendix 1 

Consider the following problem: find the most general bidyadic M e F 2 E 2 satisfying 

a A M [ck = (92) 

for all possible one-forms ex. The same problem is solved in [10] (Proposition 3.1). 
Applying the identity [18] 

aJ(o: An) = a A (ajn) + (a|a)n (93) 

valid for any one-form ex, two-form II and vector a, by choosing a so that a|a ^ 0, we 
can write from _ 

aj {ex A M [a) = ex A (aj M [a) + (a|a)M [ex = 0, (94) 

the representation 

=. = T = T a I Ml a . . 

M [ex = ex A A T , A T = (95) 

a\ot 

in terms of some dyadic A e EiFi. Since such a dyadic satisfies 

a A A T |a = M|(a A a) = 0, (96) 

we must have 

~A T \a = a\~K = Acx. (97) 

This means that a is a left eigen-one-form of the dyadic A. Because this must be valid 
for any one-form et, the dyadic is actually a multiple of the unit dyadic, 

t = AT, (98) 

and 

M[cx = AaAF = A\W T [a. (99) 
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Since M must satisfy 

(M - M^ T ) [a = (100) 
for any one-form a, by choosing basis one- forms a = e,, we finally obtain 

M - AP^ T = (M - AJ( 2 ) r ) | ^ ey ey = - Al^) \ £j e i3 = 0, 

whence M must be a multiple of the unit bidyadic . 
Multiplying (92) by e^L and applying 

e N [(uA^[u) = -v\ (e^M) \y = U m \\uu, (101) 
we obtain the following rule: if 

M m [[uv = (102) 
is satisfied for all v, the metric bidyadic must be of the form 

M m = Me w J {2)T - (103) 

Appendix 2 

Let us consider a bidyadic 

M = nC + AD + aT (2) , (104) 

defined by two two-forms II, A, two bivectors C, D and a scalar a. Let us find the inverse 
bidyadic by starting from the ansatz 

M" 1 = nC + AD' + a'T (2) , (105) 
where C, D and a are replaced by the unknown quantities C, D', a'. 
Inserting (104) and (105) in 

M|M _1 = T (2) , (106) 

and assuming that II and A are linearly independent and 

a' = 1/a, (107) 
we obtain a relation between the unknown and known bivectors as 

n|C + a A|C \fC'\ ,/C\ nn „, 

n|D A\n + a ){n> ) = ~ a [n )■ (108) 

This can be solved as 

D = (n|C + a)(A|D + a) - (A|C)(n|D). (110) 

It is easy to verify that (106) is satisfied by (105) with (109) substituted. In the case 
D = there is no inverse. 

For the special case D = in (104) the result (109) is reduced to 

C = /TT ,~ 1 -C. (Ill) 

a{n\C + a) v ; 
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Appendix 3 

Since (75) contains P-medium and Q-medium solutions, one may ask how to distinguish 
whether a given solution M m corresponds to a P-medium or a Q-medium? 

Let us assume that the modified medium bidyadic M m satisfies (75) and study first the 
Q-solution by making contraction by a one-form a as 

M m [L«« = MW 2) [[aa = M(§||aa)§- (§|a)(a|§). (112) 
Its double-cross square yields 

(M m LL««) (2) = M\^\ aa )m\<x a )W 2) -^((§|a)(a|§))), (113) 
while the double-cross cube becomes 

(M m [Lc*ao (3) = M 3 m aa ) 2 m\ aa )W 3) -W^m^m). (m) 

Applying dyadic rules [18], this can be expressed in the form 

(M m LL««) {3) = M 3 (^||«a) 2 ((^||«a)^ 3 )-^ 2 )^((^|«)(«|^))) 

= M 3 {^\\aa) 2 (e N e N A) )(e N eN[[oia). (115) 

Let us now do the same operations for the P-medium bidyadic. 

M m LL«a = M{e N ^ T )[[aa 

= Me N [(a A (P T \a) A W r ). (116) 

Denoting the simple bivector by 

e JV L(aA(P r |a))=aAb, (117) 

we have 

M m [[aa = (aAb)[P T = b(a|P T ) - a(b|P T ) (118) 
Its double-cross square yields 

(M m LL««) (2) = (b A a)((a|F) A (bp)), (119) 
and the double-cross cube vanishes 

(M m LLaa)( 3 ) = 0, (120) 

for any one-forms a. 

Now this could serve as a test for a modified medium bidyadic M m which satisfies (75) and 
is of full rank, whence both Q and P must be of rank 4. If M m satisfies (M m |_|_a:a:)( 3 ) ^ for 
any one-form ex, it corresponds to a Q-solution. In the converse case, it either corresponds 
to a P-solution or a special Q-solution with antisymmetric dyadic Q. e N e N \\Q^ = is 
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ruled out by Q being of full rank. Q antisymmetric implies also (M m |_|_ Q:Q < ' 2 ' ) = f° r a ^ 
a, but the same requirement for the P-medium bidyadic would lead to either a A b = 
or to (a|P T ) A (b|P T ) = (a A b)|P( 2 ) T = 0, which are the same condition for full-rank P. 
This condition equals o: A (P T |a) = for all a which is only valid for P being a multiple 

of T. This last case, again, corresponds to M m LL ^ = f° r a ^ a > which applied to the 
Q-medium bidyadic is contrary to the full-rank assumption. 

To conclude, we have the following test for the medium bidyadics of Q- and P-media 
satisfying (75) when the modified medium bidyadic M m is of full rank. 

• For (M m |_L Q;Q 0^ 3 ' ) 7^ we have a Q-solution 

• For (M m [[aa)^ = and (M m |_|_a") (2) ^ we have a P-solution 

• For (M m [[aa.)^ = and M m |_|_c*cx ^ we have a Q-solution with antisymmetric 
dyadic Q 

• For M m L|_a:a: = we have a P-solution with P = PL 

This result has the following conclusion. When the medium bidyadic M m is not restricted 
by a dispersion equation, i.e., it satisfies (M m [[i/uY 3 ^ = for any one-form v, it has a 
Q-medium solution only when Q is an antisymmetric dyadic, i.e., it belongs to case 1 
solutions. In other words, there are no Q-medium solutions of full rank in case 2. 
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